Boson stars in massive dilatonic gravity 
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We study equilibrium configurations of boson stars in tlie framework of a class scalar-tensor theories 
of gravity with massive gravitational scalar (dilaton). In particular we investigate the influence of 
the mass of the dilaton on the boson star structure. We find that the masses of the boson stars in 
presence of dilaton are close to those in general relativity and they are sensitive to the ratio of the 
boson mass to the dilaton mass within a typical few percent. It turns out also that the boson star 
structure is mainly sensitive to the mass term of the dilaton potential rather to the exact form of 
the potential . 
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I. INTRODUCTION 

Boson stars are gravitationally bound macroscopic 
quantum states made up of scalar bosons. They dif- 
fer from the fermionic stars in that they are only pre- 
vented from collapsing gravitationally by the Heiscn- 
berg uncertainty principle. Boson stars were first con- 
sidered by Kaup and then by Ruflini and Bonazola 
[^. They found that the boson stars as described by 
non-interacting massive, complex scalar field had masses 

of the order of w :^£pi -v^j^ere Aipi is the Planck 

niB ^ 

mass and is the boson mass. In a latter work Colpi, 
Shapiro and Wasserman ||] analyzed the consequence of 
switching on a quartic self-interaction for the boson field. 
This results in a drastic increase of the mass of the bo- 
son star, which even for small values of the coupling con- 
stant turns out to be of the order of the Chandrasekhar's 
mass when the boson mass is similar to a proton mass. 
Thus, the boson stars arise as possible candidates for 
non-baryonic dark matter in the universe. Although we 
still have no astrophysical evidences for their existence, 
the boson stars are good model to learn more about the 
nature of strong gravitational fields not only in general 
relativity but also in the other theories of gravity. 

The most natural and promising generalizations of 
Einstein's general relativity are scalar-tensor theories of 
gravity. The most studied class of such theories are 
scalar-tensor theories with massless gravitational scalar 
field including the Brans-Dicke theory as a special case. 
It should be noted that an interesting specific class of 
scalar-tensor theories with dilaton field arise from the 
low energy limit of string theory, which confirms their 
importance for current physics. 

Boson stars in the framework of scalar-tensor theo- 
ries with massless gravitational scalar have been studied 
widely, too. The first scalar-tensor model of a boson star 
was studied by Gundersen and Jensen Q| , who considered 



Brans-Dicke theory with lobd = 6. Their work was gen- 
eralized by Torres , who studied boson stars in scalar- 
tensor theories with nonconstant LUBoi^)- The conclu- 
sion is that boson stars can exist in any scalar-tensor 
theory of gravity with massless gravitational scalar, as 
the masses of the boson stars are always smaller than in 
the general relativistic case irrespective of the coupling 
functions ujbd{^)- We should mention also the interest- 
ing paper by Z. Tao and X. Xue |^ where boson stars are 
studied in the framework of a scalar-tensor theory with 
a coupling between the dilaton and the mass term of the 
boson field. 

More recently boson stars in scalar-tensor theories have 
been studied in the papers by Torres et al p, ||] and in 
the paper by Comer and Shinkai jsj. In Qboson stars 
have been studied in connection with so-called gravita- 
tional memory pO[ |, while their stability through cosmic 
history has examined in ||] and The dynamical evo- 
lution of boson stars has been investigated in the paper 
by Balakrishna and Shinkai Charged boson stars in 
scalar-tensor gravity with massless gravitational scalar 
have been studied in [Q. 

As we already mentioned boson stars in scalar-tensor 
theories have been studied only in the case of scalar- 
tensor theories of gravity with massless gravitational 
scalar. From a field-theoretical point of view it is to be 
expected that the scalar part of the gravitational field is 
massive and leads to a force with a finite range. It also 
should be noted that there is no symmetry which forbids 
a mass term, or a more general potential term for the 
gravitational scalar. Moreover, if string theory and it's 
low energy limit are relevant to the real world then the 
dilaton must be massive. Unfortunately, our current un- 
derstanding of how the dilaton acquires mass is rather 
primitive and it's tied to our lack of understanding of 
supersymmetry breaking. Since we don't have a definite 
model for dilaton mass generation the mass of the dila- 
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ton and the form of its potential are completely unknown 
from theoretical point of view. 

The purpose of the present paper is to study boson 
stars in the framework of some class of scalar-tensor the- 
ories of gravity with massive gravitational scalar. In par- 
ticular we study the possible influence of the mass and 
more generally of the potential of gravitational scalar on 
the equilibrium configurations and stability of the bo- 
son stars. The obtained results, however, are typical for 
wider classes of scalar-tensor theories with massive dila- 
ton not only for the class we consider. 



II. A SCALAR-TENSOR THEORY OF GRAVITY 
WITH MASSIVE GRAVITATIONAL SCALAR 

The most general scalar-tensor theory of gravity in- 
cluding such a term is described by the following action 



+ + Amatteri^ ■matter, 9 (1) 



where 5^,^ is the space-time metric, R is the Ricci scalar 
curvature with respect to the space-time metric. The 
functions and U ($) are, in general, arbitrary differ- 
entiable functions of the gravitational scalar <i>. A new 
method for determination of these functions using astro- 
physical observations is developed in the recent articles 
[ fi3[ , but at present we have no good enough data to 
make use of this method. The last term in denotes 
the action of the matter, which is a functional of the mat- 
ter variables, collectively denoted by ^matter, and of the 
space-time metric g^i^. 

For our purpose in this article the action (]l|) is too 
general. That's why we restrict ourselves to a subclass 
of scalar-tensor theories which don't include explicitly 
a kinetic term for the gravitational scalar in the action 
(i.e. with function = 0). In particular we consider 
scalar-tensor theories described by the action 
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d^xy/^i^^R + Ui^) 

~^ -^matter matter 1 Qi-lu 



(2) 



Such scalar-tensor theory with only one unknown func- 
tion is much more definite. It has been considered 
at first by O'Hanlon |jl4|in connection with Fujii's the- 
ory of massive dilaton It may be considered also 
as a part of more general theories of gravity like (4-|-l)- 
dimensional Kaluza-Klein model described by Fujii in 
[^6|, or model of gravity with violated local conformal 
symmetry in affine-connected spaces which probably may 
be related with string theory . A recent development 
of model of dilatonic gravity with action (^) related with 
newest astrophysical data on cosmological constant prob- 
lem may be found in pq| . In addition there it was shown 



that this model is consistent with extremely high preci- 
sion with all known data in solar system gravitational 
experiments if one chooses properly the potential U{^). 
Under proper normalization of this potential solar system 
end Earth surface experiments actually give restrictions 
only on the mass term in the potential C/($). 
The corresponding field equations are 



1 c/($). 



2 $ 



(3) 



For the treatment of the local Earth surface, astrophys- 
ical and astronomical problems in star systems (but not 
for cosmological problems in the scales of universe) we 
demand that a weak-field approximation should be pos- 
sible assuming in addition that the space-time is asymp- 
totically flat. Then in weak-field approximation we may 
write 



9 ^" ^^fiu: 



(4) 



where rj^^ is the flat space-time metric, $00 is the back- 
ground value of the $, and Shf^^, and are small per- 
turbations due to the local masses. Consistency of the 
fleld equations then requires that 



^^(1>oo) = ^(*oo) = 0. 



(5) 



Using now the consistency conditions (|5|) we obtain that 
the mass of the gravitational scalar $ is 



m 



^00 d^U 
3 



(6) 



Further on, without loss of generality we take — 1. In 
this case, taking into account that the range of the grav- 
itational scalar is flnite we obtain that the bare gravita- 
tional constant coincides with the background physi- 
cal gravitational constant. 

The PPN formalism for scalar-tensor theories of grav- 
ity with a massive gravitational scalar has been devel- 
oped by Helbig who has also calculated the basic 
gravitational effects as perihelion shift and so on. Using 
Helbig's results and the experimental data published in 
pO| and [|l| we obtain that the scalar-tensor theory we 
consider here is consistent with the experiments when the 
Compton length of the gravitational scalar $ is less than 
2 mm, or equivalently when the mass of the gravitational 
scalar satisfies the constraint m > 10~^eV |l8[ . 

So far our considerations have been made in the Jordan 
frame. For some purposes, in particular for numerical 
calculations, the scalar-tensor theories are better formu- 
lated in the conformal Einstein frame. The metric in the 
Einstein frame is given by the conformal transformation 
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(7) 



It should be stressed that it's the Jordan frame which is 
the physical frame while the Einstein frame is a conve- 
nient mathematical tool. Including also the dilaton field 
If via the formula — exp{—2aip), where a = the 
action (||) written in terms of and ip reads 



A 
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~^ •^matter matter i 



Here R is the Ricci scalar curvature with respect to 
the metric g^,v, U{(p) = A^{(p)U{^{ip)) is the poten- 
tial of the dilaton field (p and A^^ip) — ^~^{^p). It 
should be noted that the convention = 1 in terms 
of the dilaton field writes = 0. As must be ex- 
pected the dilaton potential U(}p) satisfies the conditions 
U{ipoc) — rfy^'' iVoo) — and the mass of the dilaton co- 
incides with the mass of the gravitational scalar $ which 
gives another equivalent representation of the same phys- 
ical object. 



III. BOSON STAR MODEL 



trace. The Einstein frame energy-momentum tensor T^^, 
is given by 

T^, = ^A^^) + d^^d,^+) - (12) 

^A\^) {dp^+dP^ - 2A^{^)W{^+^)) g^,. 

We consider a static and spherically symmetric boson 
star. Then the metric g^^i, can be written in the standard 
form 



We also demand a spherically-symmetric form for the 
boson field and we adopt a form consistent with the static 
metric: 



(14) 



Using the metric (|T^) and the equation defining the 
form of the boson field together with the Einstein frame 
energy-momentum tensor ( [T^ ) in the field equations ( |Tl| ) 
we get the equations for the structure of the boson star. 
Before we explicitly write them we are going to introduce 
a dimensionless radial coordinate by 



We take the matter action to be the action of complex, 
massive and self-interacting scalar field ^' which in the 
physical Jordan frame has the form 

Amatter^ jd^X^g + - W^(*+*)j (9) 



where W{^+^) 
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*+\I'-|-|;A(\I'+vE')2 is the potential 



of the boson field. Then the action of the gravity and 
matter in Einstein frame is 



A = - 
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d'^xy^ {R - 2g''''d^ipd^p + U{ip)) 



-J d'xV^(^^A'{ip)g^''d^^+d.^ 
-A'^{p)W{^+'i) 



(10) 



Varying the action ( p^ ) with respect to g^j,, p, and 
\E'+ we obtain the field equations 

G^^ = K*Tf,u + 2df,pd^p - g^^ydppdPp + ^U{p)gp^, (11) 



\/p\/Pp+-U'{p) 



VpV* + 2adppdP'^ = -2A\p>) 



2. ,dT^(«'+1') 



where Vp is the Levi-Civita connection with respect 
to the metric gpu, Tpu is the Einstein frame energy- 
momentum tensor of the boson filed ^' and T is its 



r = niBP- 



(15) 



From now on, a prime will denote a differentiation with 
respect to the dimensionless coordinate r. We also define 
other dimensionless quantities by 



= , fj — -v/k^it, A 

niB 



A 
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and dimensionless potential V{p) given by 
U{p) = m^Vip). 



m 

16 

rriB 



(17) 



With all these definitions, the equations of the structure 
of the boson star in Einstein frame reduce to the following 

A' = + re^TS + + ]--i\e^V{p^), (18) 

r 2 



1 



re li + rp 



1 



-Yre^V{p), 



2 

v' -\' 



2e^r^-2p'--i'e^V{p), 



dV{p) 



2 ' 



A' 



- I (t' - n~e 

r 



dp 



2A^{tp)e 



da^ ' 



Here 



-^T.f is the dimensionless Einstein frame 

ml J' 



energy-momentum tensor and T is its trace. In explicit 
form we have 
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T,' = \n'A\^)e-^a^ + \AH^^)e-\"+A\^)W{a^), (19) 

Ti = -\n^A^{^)er'^a^ + \A^{^)e-^a''+A\^)W{a^), 

T=-n^A^{ip)e-''a^ + A\ip)e-^a'^ + AA^{ip)W{a^). 

By reasons imposed by the numerical method used in 
the present paper, instead to investigate the system jl^ ) 
we solve numerically an equivalent system obtained as 
follows. First making use of the second equation of ( p^ ) 
we express as a function of v, v' , a, a' , if, (p' and then 
substitute in the other equations. In this way we obtain 
the system: 



z." = + + - + - YV{^)- 



the only quantity which changes under transition to Jor- 
dan frame is v. The non-singularity at the center and 
the asymptotic flatness in Jordan frame require corre- 
spondingly that v'{0) = and v{oo) = 0. Taking into 
account that v — v -\- 2aLp and Einstein frame boundary 
conditions for the dilaton field Lp we obtain that Jordan 
frame boundary conditions for v are satisfied if and only 
if = v{oo) = 0. 



IV. CONSERVED QUANTITIES 

The action of the boson field is C/(l)-invariant under 
the global gauge transformations ^' e'^^I', where a is 
a constant. This global C/(l)-symmetry gives rise to the 
following conserved current in Jordan frame 



(22) 



The same current written in terms of the Einstein 
frame metric is 



r \ r Z 4 dip 



2ap'a' + + 2A^{ip) — rV^fi- 

r \ r a<T^ 



n^e-'^a + ^iT.' + T^+^'Viip)) 



(20) 



where is given by 



I + rv' - r^ip'^ - ^A^iipya'"^ 
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-r^{\"fW{ip)-^n^A^{ip)e-''a^ + A'^{ip)W{a^)) 



(21) 



We solve numerically the system (|20|). The boundary 
conditions for the system are the following. We demand 
asymptotic flatness which means that i/(oo) = 0. The 
non-singularity of the v at the center of the star re- 
quires that J^'(O) = 0- Concerning the dilaton ip, the 
non-singularity at the center implies ip'{0) = 0, while 
at infinity we must have (p{oo) = as required by the 
asymptotic flatness. Non-singularity of the boson fleld 
at the center implies a'{0) — and the asymptotic flat- 
ness requires that a{oo) = 0. To complete the nonlinear 
eigenvalue problem for we also must give the central 
value of the boson fleld ct(0) = Cc- 

Here we must make some comments. We have im- 
posed the boundary conditions in Einstein frame. But 
what is important is the non-singularity and asymptotic 
flatness in the physical Jordan frame. However, it's not 
difficult to see that all Einstein frame boundary condi- 
tions phrased in Jordan frame remain the same. Indeed, 



.r - -A2(^)V^5^'' {^>d,^+ - VI/+9,*) . (23) 

The conserved current leads to conserved charge - the 
total number of particles making up the star: 



N : 



(24) 



It should be stressed that N is the total number of par- 
ticles in the physical Jordan frame. Taking into account 
the explicit form of the metric and of the boson field we 
obtain 



(25) 



In order to calculate numerically the total number of 
particles we must present the above integral in terms of 
Einstein frame metric. This can be done easily by pre- 
senting the Jordan frame metric as conformally trans- 
formed Einstein frame metric. This way we obtain 

Mi, = msN = ( Q H drr^A\^)e^a^ (26) 
V 2mB ) Jo 

where Aifi \s the rest mass of the star in the Jordan 
frame. 

The binding energy in Jordan frame is then defined by 



£b = M-Mi 



(27) 



where M. is the total mass of the star. 

In contrast to the general relativity the definition of 
mass in scalar-tensor theories is subtle . There are three 
possible mass definitions in Jordan frame, notably the 
Schwarzschild mass Ms (i-e. ADM mass in Jordan 
frame), the Kepler mass M.k and the tensor mass Mt 
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which is the ADM mass in Einstein frame [g2|. These 
mass definitions are equivalent in general relativity, but 
they may give difi'erent results in scalar-tensor theories 
of general type. It's the tensor mass which has got phys- 
ically acceptable properties: The tensor mass has the 
important property to peak, as a function of the central 
density at the same location where the particle number 
(respectively the rest mass) takes its maximum. This 
property has been proved in the most general scalar- 
tensor theory including a potential term for the gravi- 
tational scalar (dilaton), see for details ||2^. That the 
tensor mass and particle number peak at the same lo- 
cation results in a cusp in the bifurcation diagram M 
(respectively £b) versus Mr. 

All masses we have mentioned are asymptotically mea- 
sured quantities. Taking into account that in the scalar- 
tensor theory we consider the gravitational scalar (dila- 
ton) varies essentially in a finite space-time domain , it's 
not difficult to see that all masses in our case are identi- 
cal. That's why all properties proven for the tensor mass 
are shared also by the others and the mass of the star may 
be expressed in three different ways as one wishes. We 
take the expression for the mass of the star in the form 
of the Schwarzschild mass in Einstein frame, namely 



M 



For numerical purposes we also introduce the following 
dimensionlcss masses 



M = M/ 



PI 



M2 



2mB ^niB 
Therefore the dimensionlcss binding energy is 

Eb = M- Mr . 



(29) 



(30) 



We also define an effective dimensionlcss radius of the 
boson star as: 



R 



N 



(31) 



V. NUMERICAL RESULTS AND DISCUSSIONS 

To solve the nonlinear ODE's (|o|) with the correspond- 
ing boundary conditions and spectral parameter Q we 
apply the continuous analog of the Newton method , 
^5| . In our case this method leads to a linear systems 
ODE's with respect to the increments of unknown func- 
tions coupling with an algebraic equation for the incre- 
ment of fl. The appropriate linear boundary value prob- 
lems we solve numerically using the spline collocation at 
the Gaussian points on the irregular mesh condensing to 
the centre r — of the star. For details see our work pq . 

As we have already mentioned at present we have only 
some weak experimental constrains on the mass of the 
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FIG. 1. The dependence of the mass and rest mass of the 
boson star on the central value ac 



dilaton and on the form of its potential. A definite theo- 
retical suggestions for them do not exist. That's why we 
consider first the simplest form for the potential of the 
gravitational scalar 



-m2($-l)2. 
2 ^ ' 



[/($) 

The corresponding dilaton potential is 

U{^) = mV(^) = ^m' (^1 - exp(^) 



(32) 



(33) 



The numerical results for the chosen potential are pre- 
sented in the figures. 

Fig.l presents a configuration diagram which is typical 
for large variations of the parameters not only for those 
values presented on the figure. 

The dependence of the boson star mass on the central 
value cTc is presented in Fig. 2 for three different values 
of the ratio 7 = at fixed value of the parameter 
A — 10. It's seen that the curves M{(Jc) differ slightly 
from each other. More detailed investigation for much 
more values of 7 shows that the curves tend to the same 
curve M{ac) from the general relativistic case for large 
values of 7 (7 « 10). This behaviour is confirmed also 
by Fig. 3 where the dependence of the boson star mass on 
the parameter 7 is shown for Gc = 0.1 and A = 10. For 
small values of 7 the boson star mass increases when the 
parameter 7 increases. Then for larger values of 7 the 
mass M tends to a fixed value which is just the general 
relativistic value. The cause for this behaviour is that 
the increasing of the parameter 7 is actually an increas- 
ing of the dilaton mass with respect to the boson mass 
which leads to smaller dilaton range and therefore to a 
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FIG. 3. The dependence of the boson star mass on the 
parameter 7 for A = 10 



The dependence i/(r) is shown in Fig. 5 for three dif- 
ferent values of 7. As it is seen the influence of 7 on i/{r) 
is slight. 




suppressing of the dilaton field which means that the gen- 
eral relativity is recovered. The Fig. 3 also shows that 
the boson stars in the model under consideration always 
have masses smaller than those in general relativity. The 
deviations from general relativity are within a typical few 
percent. 

The influence of the dilaton mass, respectively 7, on 
the dependence of the boson field on the radial coordi- 
nate r is very slightly as one may see from Fig. 4. The 
dependence a{r) is presented there for three different val- 
ues of 7. The three curves are extremely close to each 
other. 



FIG. 5. The physical function ^(r) for three different val- 
ues of the parameter 7. 

The most sensitive to the value of the parameter 7 is 
the dilaton field cp. The dependence (p{r) is presented 
in Fig. 6 for three different values of 7. One sees that 
the increasing of 7 leads to suppressing the dilaton field 
as one must expect. The dependence 3>(r) is shown in 
Fig. 7. At the center of the star one has <&(0) > 1 
which means that the gravitational constant at the center 
satisfies G(0) < 1. Therefore the boson stars are less 
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FIG. 6. The dependence ip{r) for three different values of 
the parameter 7. 



sensitive mainly to the second derivative of the dilaton 
potential ^^-^^t^{0) determining the dilaton mass rather 
to the exact form of the potential. 
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FIG. 7. The function $(r) for three different values of the 
parameter 7 

gravitationally bounded than in general relativistic case. 

We have also examined the cases of different values of 
A. The picture is qualitatively the same. 

It is interesting to know how the exact form of the dila- 
ton potential V{ip) influences the boson star structure. 
We have examined several potentials V{ip) with the same 
dilaton mass (for example V{(p) = 2(p^, V{ip) = 2sm^{ip), 
V{ip) = 4(l-cos((^)), V{ip) = 2{l-cxp{-Lp^))). Our nu- 
merical results show that boson star structure in practice 
does not depend on the exact form of the dilaton poten- 
tial in the class of potential we have examined. This 
gives strong evidences that the boson star structure is 



FIG. 8. The boson star mass Ad as a function of the boson 
star radius R 

Finally in Fig. 8 one can see that the usual dependence 
of the mass M{R) of boson star on its effective radius R 
take place in the scalar-tensor model with massive dila- 
ton, too. As seen, there exist a domain of stability of the 
boson star, followed by domains of unstable states. 

VI. BOSON STAR STABILITY 

Here we briefly discuss stability property of the boson 
stars we consider. Our analysis is based on catastrophe 
theory 

The basis for our analysis is Fig. 9 where the binding 
energy Ei, is plotted against the rest mass Mn (i.e. parti- 
cle number). This figure, as our numerical results show, 
is typical for a large number of values of 7 and A. There- 
fore the forthcoming conclusions concern the general case 
not only the special case presented by the figure. 

Fig. 9 is actually a bifurcation diagram. For small cen- 
tral densities the binding energy is negative which shows 
that the boson star is potentially stable and one assumes 
that the boson star is stable against small radial pertur- 
bations. As the central density is increased one meets 
a cusp. The second branch as a whole has higher mass 
and therefore it is unstable. At the cusp, the boson star 
stability changes - one radial perturbation mode develops 
instability. 
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FIG. 9. The boson star binding energy Eb as a function 
of the rest mass Mr 



boson stars, our results are qualitatively close to the re- 
sults obtained for boson stars in Brans-Dicke theory 

, and for boson stars considered in a scalar-tensor the- 
ory with a kinetic term and simple coupling between the 
dilaton and the mass term of the boson filed in Jordan 
frame |^ which in our notations reads TO^<i>^+^'. In 
the domain of instability of boson stars the results of the 
last reference differ essentially from ours. As in general 
relativity and Brans-Dicke theory beyond the point of 
stability for large central densities the boson star mass 
in our model of massive dilatonic gravity drops, oscillates 
a bit and approaches a constant value while in the model 
considered in ||^ it increases rapidly. 
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VII. CONCLUSION 



In this paper we have analyzed static spherically sym- 
metric boson stars in a scalar tensor theory of gravity 
with massive gravitational scalar (dilaton). We have 
studied their equilibrium properties for different values 
of the ratio 7 = -^2- . The conclusion is that stable boson 
stars may exist for any value of 7. When 7 increases the 
equilibrium configurations tends to those from general 
relativity and for sufficiently large 7 the general relativis- 
tic case is in practice recovered. The masses of the boson 
stars are turned out to be always smaller than those in 
general relativity. The typical deviations of the masses in 
our model of dilatonic gravity from those in general rel- 
ativity is a few percent. The small deviations of bosonic 
star structure in the considered model of dilatonic grav- 
ity from general relativistic case are significantly greater 
then the corresponding deviations in solar system and 
Earth surface experiments studied in It turns out 

also that the boson star structure is sensitive mainly to 
the mass term of the dilaton potential rather to the exact 
form of the potential. The boson stars in our case are less 
gravitationally bound as a consequence of the fact that 
the gravitational constant G within the stars is smaller 
than one. 

For completeness we have examined also the case of 
scalar-tensor theories with a massive gravitational scalar 
with kinetic term ujbd^~^ g^'^d^^^d^^. In these cases 
the obtained results for different values of the parameter 
liJbd are qualitatively the same as in the case ojbd = Q 
and for large lobd the differences between bosonic stars 
in dilatonic gravity and in general relativity are even 
smaller. 

It is worth noting that, in the domain of stability of 
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